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Abstract. While the modern stellar IMF shows a rapid decline with increasing mass, 
theoretical investigations suggest that very massive stars (> 100 Mq) may have been 
abundant in the early universe. Other calculations also indicate that, lacking metals, 
these same stars reach their late evolutionary stages without appreciable mass loss. 
After central helium burning, they encounter the electron-positron pair instability, 
collapse, and burn oxygen and silicon explosively. If sufficient energy is released by 
the burning, these stars explode as brilliant supernovae with energies up to 100 times 
that of an ordinary core collapse supernova. They also eject up to 50 Mq of radioactive 
56 Ni. Stars less massive than 140 Mq or more massive than 260 Mq should collapse 
into black holes instead of exploding, thus bounding the pair-creation supernovae with 
regions of stellar mass that are nucleosynthetically sterile. Pair-instability supernovae 
might be detectable in the near infrared out to redshifts of 20 or more and their ashes 
should leave a distinctive nucleosynthetic pattern. 


1 Introduction 


Owing to the lack of any metals, the cooling processes that govern star formation 
are greatly reduced for first generation of stars (Pop III). Magnetic fields and 
turbulence may also be less important at these early times jl|. Consequently, 
theoretical studies (TTJ indicate that the Jeans mass for primordial stars in their 
special environment may have been as great as ~ 1000 Mq. Numerical simulation 
of primordial star formation predict the occurrence of such stars at red shifts 
~ 20 and an initial mass function (IMF) that either peaks at ~ 100 M 0 jj],[]J or 
is bimodal |l^|, i.e., also contains stars of a few Mq. 

Once formed, at solar metallicity, massive stars ordinarily experience signifi¬ 
cant mass loss || and may end as relatively small objects, but for low metallicity 
mass loss is suppressed. In § || we discuss the peculiarities of mass loss and evo¬ 
lution of very massive primordial stars. Figure [I] gives an overview of expected 
final fates of metal-free stars as a function of initial mass. In § |3| we examine 
the expected light curve of a pair-creation supernova from a 250 M 0 star at a 
redsliift of z = 20 and in § i we review nucleosynthetic yields from Pop III. Some 
conclusions are given in § 5. 
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Fig. 1. Initial-final mass function of non-rotating Pop III stars. The x-axis gives the 
initial mass. The y-axis gives both the final mass of the collapsed remnant (thick black 
curve) and the mass of the star when the event begins that produces that remnant 
(e.g., mass loss in AGB stars, supernova explosion for those stars that make a neutron 
star, etc.; thick gray curve). We distinguish four regimes of initial mass: low mass 
stars below ~ 10 Mq that end as white dwarfs; massive stars between ~ 10 Mq and 
~ 100 Mq that form an iron core that eventually collapses; very massive stars between 
~ 100 Mq and ~ 1000 Mq that encounter the pair instability; and supermassive stars 
(arbitrarily) above ~ 1000 Mq . Since no mass loss is expected for Z = 0 stars before 
the final stage, the grey curve is approximately the same as the line of no mass loss 
(dotted). Exceptions are ~ 100 — 140 Mq where the pulsational pair instability ejects 
the outer layers of the star before it collapses, and above ~ 500 Mq where pulsational 
instabilities in red supergiants may lead to significant mass loss |^j . Since the magnitude 
of the latter is uncertain, lines are drawn dashed. For a more detailed description, please 
refer to [|t|] 
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Fig. 2. Growth rate, K, as a function of the age of the star. The amplitude of the 
pulsation grows with time as A(t) = Aq exp{2^if} exp{A't} Jsj]. 


2 Mass Loss 

Current studies of winds driven from hot stars by interactions with atomic lines 
indicate that their mass loss rate decreases with metallicity as Z 1 / 2 |||| or even 
more jl6|. This scaling seems to hold down to 0.1 % solar metallicity. Extrapolat¬ 
ing to zero, it seems reasonable that mass loss from this wind driving mechanism 
becomes negligible. Winds driven by continuum opacity in low metal stars are 
not very well understood (Kudritzki, priv. com.), and will be neglected here. 

Stars above ~ 60 M 0 may also be unstable to the epsilon mechanism - pul¬ 
sations driven by the high temperature sensitivity of nuclear burning ||~i7]]T7i| . 
However, historical studies of stellar stability focused on stars of solar compo¬ 
sition and did not take into account the peculiarities of the first generation of 
stars. We have recently studied the pulsational instability due to epsilon mecha¬ 
nism in massive metal-free stars between 120 and 1000 M 0 ||. Their structure is 
uniquely different from later stellar generations. Since hydrogen burning through 
the “pp-chain” is not sufficient, the star contracts to a high enough tempera¬ 
ture to produce carbon by the triple-alpha reaction, more than 1x10 8 K in the 
stars investigated here. A mass fraction of ~lxl0~ 9 is sufficient to stop the con¬ 
traction and supply energy by the CNO cycle for hydrogen burning. The stars, 
however, stay very compact and hot in their centers throughout their hydrogen 
burning lifetime. The higher temperature causes a lower temperature sensitivity 
of the nuclear energy generation. From pulsational analysis we find that the ep- 
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silon mechanism is weak in these compact stars and operates only for a fraction 
of the hydrogen-burning life-time (Fig. |) 0 - Radial pulsations driven by opac¬ 
ity or recombination are not found. We estimate that the resulting mass loss due 
to the epsilon mechanism should be less than 5 % for a 500 M 0 star, and perhaps 
10% for a 1000 Mq star, but quite negligible for stars of lower mass. Stars of 
~ 500 Mq or more may encounter a red supergiant phase towards the end of 
central helium burning. This could result in an additional, maybe significant, 
mass loss, but its strength is not yet known. Therefore stars of < 500 M 0 can 
be safely assumed to reach carbon burning without significant mass loss. 


3 Supernovae at the Edge of the Universe 


Stars that reach carbon burning with a helium core mass of ~ 64... 133 M 0 
(corresponding to ZAMS masses of ~ 140 ... 260 M 0 for stars without mass loss) 
are unstable to pair-creation (see Fig. [j]), collapse, then explode as a supernova 
(SN). These are the most powerful thermonuclear explosions in the universe. 
Pair-creation SNe release energies ranging from 3xl0 51 erg for a 64 M 0 He core 
up to almost 100xlO 51 erg for a 133 M 0 He core 0 enough energy to disrupt a 
small proto-galaxy. In Fig. ^|we show, in the observer frame, the early light curve 
of an exploding 250 M 0 star, neglecting intergalactic and interstellar absorption. 
This star has a He core of 120 M 0 and a total explosion energy of 65xl0 51 erg, 
producing 21 M 0 of 56 Ni. The visible brightness is not quite as impressive since 
most of the energy is kinetic, but nevertheless should be a few times brighter 
than a typical Type la SN. The bolometric luminosity of an event at z = 20 is 
not much dimmer than at a red shift of a few, where Type la SNe have already 
been observed. The main difference, however, is the significantly larger red shift 
and time dilation. Given the bigger intrinsic time scale associated with the large 
ejected mass, they also last much longer. 

For a current standard cosmology (Oa = 0.7, h? m = 0.3, H 0 = 65km/s/Mpc, 
F?b = 0.02//i 2 = 0.047) and assuming that /i st = 10 -6 of all baryons goes into 
stars of Mi s t = 250 M 0 , at a red shift of 2 = 20, we estimate the pair-creation 
supernova rate by 


P2SN = 477 



7-7- PcW + ~) 3 

1 + z 


f 1st 
Mi s t 


47rrf lum C /°c^b/lst 

Mist 


where p c = 3Hq 2 /8nG, and the luminosity distance di um = 243 Gpc. This gives 
~0.16 events per second per universe, i.e., ~3.9xl0 -6 events per second per 
square degree. The first peak of the light curve (Fig. ||) lasts for about a month, 
the second peak (“plateau”; not depicted here), probably brighter in the far 
infrared, would peak for about 10 yr (based on preliminary calculations by Phil 
Pinto, priv. com.). Statistically, at any time per square degree about a dozen of 
these supernovae should be at the peak of the light curve, and more than ~1000 
in the plateau phase of the light curve. Note that di um oc z for high z and thus 
this rate depends critically on the red shift adopted as well as on the baryon 
fraction assumed to go into these stars. 
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Fig. 3. Preliminary light curve of pair-creation supernova from a 250 Mq star at 2 = 20 
as computed by the KEPLER code (Tt). Time, wave lengths and magnitudes (without 
internal or intergalactic extinction) are given in observer rest frame. Wave lengths that 
are beyond the IGM Ly-a absorption (2.55/rm) are displayed as dotted lines. “PRIME” 
and “NGST” corresponds to 3.5 and 5.0/im. The “spherically symmetric” emission has 
been folded to account for the extent of the “photosphere”. The first “bump” at ~ 10 3 s 
is from the shock breakout, the right “peak” is the peak of the SN light curve. 


4 Nucleosynthesis 

In massive stars most of neutrons responsible for the “weak component” of the 
s-process are made from initial CNO seeds during helium burning. Since Pop 
III stars are devoid of these seeds, both the s-process and the production of a 
“neutron excess” for the advanced burning stages is strongly suppressed 0 . In 
massive stars the odd -Z elements are therefore underproduced with respect to 
the a elements [[fj . Above ~ 25 Mq significant fallback onto the remnant already 
occurs after the SN Q (Fig. |l|), and few or no heavy elements can be ejected. 
From ~ 40 to ~ 100 Mq a black hole is formed directly |i) and essentially 
the whole star will be swallowed — no nucleosynthesis products are ejected. 
Between ~ 100 and ~ 140 M 0 the pulses of pair instability || will eject the 
outer layers of the star, possibly including parts of the CO core, especially at 
the high-mass end of the regime, but nothing heavier than magnesium leaves 
the star. After the pulses these objects probably encounter the same fate as 
their lighter cousins, i.e., the remaining part of the star falls into a black hole. 
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Fig. 4. Elemental production factors relative to solar in Pop III pair-creation super¬ 
novae (helium core mass 65-130 Mg). The thin and thick ends of the triangle and the 
dots correspond to IMF slopes of dlog IV/dlogM = —1.5, —4.5, and —2.5, respectively. 


For initial masses of ~ 140.. .260Mq, the pair instability completely disrupts 
the star when explosive burning of oxygen and silicon release enough energy to 
reverse the collapse into an explosion. However, this explosion is too rapid and 
the star not dense enough at the point when the implosion turns around to lead 
to significant neutronization. A very marked odd-even pattern results (Fig. Q) 
with silicon being the biggest overproduction and elements above germanium 
are essentially not produced. In stars more massive than ~ 260 M 0 , photo¬ 
disintegration of nuclei becomes important and the collapse is not reversed. The 
whole star falls into a black hole (Fig. [j]). 

5 Conclusions 

Due to the absence of metals the first generation of star will likely not experience 
significant mass loss by radiation-driven stellar winds or opacity-driven pulsa¬ 
tions. Their unique structure also prevents significant mass loss by the epsilon 
mechanism. Therefore very massive Pop III single stars reach carbon burning 
with enough mass to encounter the pair instability. Since current theoretical 
studies indicate that such star may constitute a significant, if not dominant, 
fraction of Pop III, we predict that their nucleosynthetic yields may have a 
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unique imprint on the chemical evolution of the early universe. Their produc¬ 
tion of odd -Z elements is by ~ 2 orders of magnitude lower than that of even -Z. 
Elements heavier than zinc are not produced. 

Since stars of lower mass (< 140 M 0 ) or higher mass (> 260 M 0 ) collapse 
into black holes without significant heavy element creation. Pair-creation SNe 
are thus a “clean” source of nucleosynthesis in the sense that neighboring mass 
ranges do not “pollute” the sample. In case of a bimodal Pop III IMF, massive 
stars in the range ~ 8 ... 40 M 0 will also contribute to the resulting abundance 
pattern, though on a slightly slower time-scale (factor ~ 2). Also in them, the 
lack of initial CNO “seeds” will lead to an elemental odd-even pattern, though 
much less expressed and they possibly contribute r-process isotopes. It is 
even conceivable that Pop III AGB stars contribute some s-process [[jl|. The 
interaction of the ejecta with the surrounding matter, a possible enrichment of 
the intergalactic medium and mixing of contributions from different-mass sources 
before the formation of the first Pop II stars will have to be studied in more detail 
in the future. 

We predict that Pop III pair-creation SNe might be detectable by future 
near infrared space experiments — all the way out to the edge of universe — 
to redshifts of 20 or more. Combined with the challenge to find old Pop II stars 
that show the predicted abundance pattern from the ashes of these explosions, 
this should allow deeper insight into the happenings at the times when the first 
sparks of stellar light terminated the “dark ages”. 
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